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ABSTRACT 

We combine new Wide Field Camera 3 IR Channel (WFC3/IR) F160W (H m ) imaging data for NGC 1399, 
the central galaxy in the Fornax cluster, with archival F475W (# 475 ), F606W (V 606 ), F814W (7 gl4 ), and F850LP 
(■Z85o) optical data from the Advanced Camera for Surveys (ACS). The purely optical #475-/814, V606 - 78i4, 
and g 4 75-Z850 colors of NGC 1399's rich globular cluster (GC) system exhibit clear bimodality, at least for 
magnitudes /si4 > 21.5. The optical-IR hu-H\a) color distribution appears unimodal, and this impression 
is confirmed by mixture modeling analysis. The V606 _ ^i60 colors show marginal evidence for bimodality, 
consistent with bimodality in V(,Q6-k\A and unimodality in hu-H\oa. If bimodality is imposed for /814-//160 
with a double Gaussian model, the preferred blue/red split differs from that for optical colors; these "differing 
bimodalities" mean that the optical and optical-IR colors cannot both be linearly proportional to metallicity. 
Consistent with the differing color distributions, the dependence of hu—Hiaa on #475-/814 for the matched GC 
sample is significantly nonlinear, with an inflection point near the trough in the #475-/814 color distribution; the 
result is similar for the /814-//160 dependence on #475— Zgso colors taken from the ACS Fornax Cluster Survey. 
These #475-2850 colors have been calibrated empirically against metallicity; applying this calibration yields 
a continuous, skewed, but single-peaked metallicity distribution. Taken together, these results indicate that 
nonlinear color-metallicity relations play an important role in shaping the observed bimodal distributions of 
optical colors in extragalactic GC systems. 

Subject headings: galaxies: elliptical and lenticular, cD 
clusters — globular clusters: general 



galaxies: individual (NGC 1399) — galaxies: star 



1. INTRODUCTION 

The major star formation episodes in the history of any 
large galaxy will be imprinted in the properties of the galaxy's 
star cluster population. Interpreting the observed properties to 
derive the formation histories has proven to be a difficult task. 
All giant elliptical galaxies contain large populations of glob- 
ular clusters (GCs), often numbering in the thousands (e.g., 
Harris 1991). The GC systems generally follow bimodal dis- 
tributions in optical colors (Zepf & Ashman 1993; Gebhardt 
& Kissler-Patig 1999). For giant ellipticals in clusters, the 
two peaks in the color distributions are roughly equal in size 
(e.g., Peng et al. 2006; Harris et al. 2006), except at large 
radii where the color distribution is more strongly weighted 
towards the blue (e.g., Dirsch et al. 2003; Harris 2009). 

If the optical colors are interpreted as a direct proxy for 
metallicity, then the bimodality represents an extraordinary 
constraint on the star formation histories of giant ellipti- 
cals. In this case, the two color peaks represent two distinct 
cluster populations, differing by roughly a factor of ten in 
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mean metallicity. The bulk of the host galaxy's stellar mass 
would then likely originate from two distinct major forma- 
tion episodes. Ashman & Zepf (1992) originally predicted 
such GC bimodality based on the idea that ellipticals form 
from gas-rich major mergers of late-type galaxies that already 
possessed extensive metal-poor GC populations. Because el- 
lipticals are believed to have formed in a more stochastic, 
hierarchical fashion, a number of other merger or accretion 
scenarios were later proposed to account for the observed bi- 
modality (Forbes et al. 1997; Kissler-Patig et al. 1998b; Cote 
et al. 1998; Beasley et al. 2002; Kravtsov & Gnedin 2005). 
However, these different scenarios are often difficult to distin- 
guish observationally from one another, and none appears to 
account naturally for all the data (see Peng et al. 2006). 

More recently, the assumption of optical colors as a sim- 
ple, linear proxy for metallicity has been reexamined (Richtler 
2006; Yoon et al. 2006). It has been known for many years 
that the slope of the metallicity as a function of optical color 
becomes shallower (i.e., color becomes more sensitive to 
metallicity) at intermediate metallicities (Kissler-Patig et al. 
1998a). Yoon et al. (2006) showed that the "wavy" nonlinear 
color-metallicity relation predicted by their stellar population 
models matched, at least qualitatively, the color-metallicity 
data assembled by Peng et al. (2006). They further pointed 
out that the "projection" from metallicity to color with such 
wavy relations can produce bimodal color distributions from 
unimodal distributions in metallicity. Cantiello & Blakeslee 
(2007) confirmed the Yoon et al. (2006) result, in the sense 
that other sets of models that include realistic prescriptions for 
the behavior of the horizontal branch as a function of metallic- 
ity also give nonlinear color-metallicity relations that can pro- 
duce bimodal color distributions. More recently, Yoon et al. 
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(201 la,b) find that the GC metallicity distributions derived for 
several galaxies from optical colors using their model color- 
metallicity distributions are not bimodal, but are similar to 
the metallicity distributions found for stellar halos in ellipti- 
cal galaxies. 

This issue remains highly controversial, but it is clear that 
additional studies of GC metallicity distributions are needed. 
In particular, it is important to constrain independently the 
form of the metallicity distributions in galaxies with promi- 
nently bimodal GC color distributions. Large samples of 
spectroscopic metallicities in extragalactic GC systems (more 
than a few percent of the total population) are now becoming 
available (Beasley et al. 2008; Foster et al. 2010; Alves-Brito 
et al. 2011; Caldwell et al. 2011; see Section |7]for a discus- 
sion of these results). However, at the distances of the Virgo 
and Fornax clusters, spectroscopic data of sufficient quality 
remains observationally expensive, requiring multiple nights 
of ten-meter class telescope time to amass significant samples. 

Another way to constrain GC metallicity distributions is 
through a combination of near-infrared (near-IR) and opti- 
cal photometry (e.g., Kissler-Patig et al. 2002; Beasley et al. 
2002; Puzia et al. 2002; Kundu & Zepf 2007; Kotulla et al. 
2008; Chies-Santos et al. 2011a,b). Hybrid optical-IR colors 
such as I-H or I-K for old stellar systems such as GCs are 
mainly determined by the mean temperature of stars on the red 
giant branch, which in turn depends almost entirely on metal- 
licity (e.g., Bergbusch & VandenBerg 2001; Yi et al. 2001; 
Dotter et al. 2007). Optical-IR colors that involve bluer pass- 
bands, such as V-H or B—K, also have the strong metallicity 
dependence from the giant branch temperature, but they have 
additional sensitivity to the main sequence turnoff, which de- 
pends strongly on age, and to the horizontal branch morphol- 
ogy, which behaves nonlinearly with metallicity and also de- 
pends on age (e.g., Lee et al. 1994; Sarajedini et al. 1997; 
Dotter et al. 2010). Although the dependence of giant branch 
temperature (and thus of I-H and similar colors) on metallic- 
ity is not necessarily linear, there is no evidence for a sharply 
nonlinear transition with metallicity, as occurs for the hori- 
zontal branch. 

The extensive GC system of NGC 1399, the dominant el- 
liptical galaxy in the Fornax cluster, has been a frequent tar- 
get for optical photometric and spectroscopic studies since the 
pioneering work of Hanes & Harris (1986). At a distance of 
20 Mpc (Blakeslee et al. 2009) Fornax is the second nearest 
galaxy cluster after Virgo, and NGC 1399 is at its dynami- 
cal center (Drinkwater et al. 2001). NGC 1399 was one of the 
first external galaxies reported as having a bimodal metallicity 
distribution, based on optical colors (Ostrov et al. 1993). In- 
terestingly, this galaxy provided the first indication of the im- 
portance of the detailed shape of the color-metallicity relation: 
Kissler-Patig et al. (1998a) measured spectroscopic metallic- 
ities for a sample of GCs in NGC 1399 and found that the 
slope of metallicity versus (V-I) color was significantly flat- 
ter than the extrapolation from low-metallicity Galactic GCs. 
NGC 1399 also has the largest sample of measured GC radial 
velocities to date (Schuberth et al. 2010). 

Here we present new near-IR and optical photometry of 
GCs in NGC 1399 using the Infrared Channel of the Wide 
Field Camera 3 (WFC3/IR) and Wide Field Channel of the 
Advanced Camera for Surveys (ACS/WFC) on board the 
Hubble Space Telescope (HST). We take an empirical ap- 
proach in this work, comparing the color distributions and 
examining the color-color relations, without trying to judge 
between different sets of models for the present. The follow- 



ing section summarizes the observational details and image 
reductions for both the WFC3/IR near-IR and ACS/WFC op- 
tical data. Section|3]describes our photometric measurements 
and selection of GC candidates. Section |4] presents the GC 
color-magnitude diagrams (CMDs) and color distributions re- 
sulting from our new photometry, as well as mixture model- 
ing analysis results for the different distributions. The form 
of the color-color relation between I%\h-H\(,q and #475-/814 is 
discussed in detail in Section|5] In Section[6] we cross-match 
our photometry against the ACS Fornax Cluster Survey (here- 
after, ACSFCS; Jordan et al. 2007) catalogue for NGC 1399 
and discuss the implications for the underlying metallicity dis- 
tribution. The final section places our study in the larger con- 
text of GC color and metallicity studies, before listing our 
main conclusions. 

2. OBSERVATIONAL DATA SETS 

NGC 1399 was observed for one orbit in the F110W and 
F160W bandpasses of WFC3/IR in 2009 December as part 
of HST program GO-1 1712. The WFC3 IR Channel uses a 
1024 2 pix HgCdTe detector with an active field of view of 
<~ 2.'05 x 2.'27, and a mean pixel scale of about 0"128 pix" 1 
(see Dressel et al. 2010 for more information). A primary 
goal of GO-1 1712 is to obtain an empirical calibration for 
the surface brightness fluctuations (SBF) method in these two 
WFC3/IR bandpasses similar to those derived for ACS in 
F814W and F850LP (Blakeslee et al. 2009, 2010b; Mei et al. 
2005). To this end, GO-1 1712 targeted 16 early-type galaxies 
in the Fornax and Virgo clusters over a wide range in luminos- 
ity and color. Another important goal of the project is to in- 
vestigate the optical-IR colors of the GC populations in these 
galaxies, using primarily the F160W bandpass (which affords 
the wider baseline) combined with existing ACS optical data. 
NGC 1399 was one of the first large galaxies observed in this 
project, and we have used it to optimize our image process- 
ing pipeline. Analyses of the SBF properties and GC colors 
for the full sample of GO-1 1712 galaxies will be presented in 
forthcoming works. 

A total of 1 197 s of integration was acquired with WFC3/IR 
in the F160W band. The data were retrieved multiple times 
from the STScI archive as the calibration reference files for 
WFC3 were updated. The photometric results improved 
markedly after flat fields produced on-orbit (described by 
Pirzkal et al. 201 1) were implemented in the STScI pipeline. 
The results presented here are derived from images retrieved 
from STScI in 2011 January. We combined the individual 
calibrated WFC3/IR exposures into a final geometrically cor- 
rected image using the Multidrizzle (Koekemoer et al. 2002) 
task in the STSDAS packag^J After experimenting with dif- 
ferent Drizzle (Fruchter & Hook 2002) parameters, we settled 
on the square interpolation kernel with a "pixfrac" value of 0.8 
and a final pixel scale of 0."1 pix" 1 . This scale is convenient 
because it is exactly twice that of ACS/WFC. 

NGC 1399 has been observed many times before with 
HST, including several times with the ACS/WFC, allowing 
us to investigate the optical-IR colors of its GCs. Cali- 
brated observations in F606W from GO- 10129 (PI: Puzia) 
and F475W+F814W from GO-10911 (PI: Blakeslee) were 
retrieved from the STScI archive and processed with Apsis 
(Blakeslee et al. 2003) to produce summed, geometrically cor- 

8 STSDAS is a product of the Space Telescope Science Institute, which is 
operated by AURA for NASA. 
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Fig. 1 . — HST WFC3/IR F160W image of NGC 1399 (top), and the same following galaxy isophotal model and sky subtraction (bottom). The 
field size is approximately 2'.3 x 2'.l, and the image is shown at the observed orientation. Pixels flagged as bad in the data quality arrays are 
zeroed in both panels. The central l"2 of the galaxy was not modeled, and is set to zero in the lower panel. 
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rected, cosmic-ray-cleaned images for each bandpass. Reduc- 
tion of the GO- 1091 1 imaging data is described in more detail 
by Blakeslee et al. (2010b); the same procedures were applied 
for the GO- 10129 program, which observed nine contiguous 
fields in F606W. We processed only the four pointings that 
overlapped with the ACS/WFC images in other bands; three 
of these pointings overlap with our smaller WFC3/IR field. 
NGC 1399 was also observed in F475W and F850LP as part 
of HST program GO-10217 (Jordan et al. 2007); in this case, 
we used the photometry catalogue produced by that program, 
as discussed in Section[6] 

Throughout this study, we employ the natural photomet- 
ric systems defined by the instrument bandpasses, rather 
than converting to the Johnson system. We calibrated the 
WFC3/IR photometry using the AB zero-point coefficients 
given by Kalirai et al. (2009) and the ACS photometry us- 
ing the AB zero points from Sirianni et al. (2005). We cor- 
rected for Galactic extinction towards NGC 1399 assuming 
E(B-V) = 0.0125 mag (Schlegel et al. 1998), the ACS/WFC 
extinction ratios from Sirianni et al. (2005), and the //-band 
extinction ratio from Schlegel et al. (1998). Table Q] sum- 
marizes the observational details of the data sets used in the 
present study, including in the last column the symbols used 
to denote magnitudes in the various bandpasses. 

3. OBJECT SELECTION AND PHOTOMETRY 

In order to obtain photometric catalogues of GC candidates, 
we first constructed elliptical isophotal models for the differ- 
ent bandpasses as described in our previous works (e.g., Tonry 
et al. 1997; Jordan et al. 2004; Blakeslee et al. 2009, 2010b) 
and used these to subtract the galaxy light. Figure Q] shows 
our WFC3/IR F160W image of NGC 1399 before and after 
the galaxy subtraction. The compact sources corresponding 
to GC candidates, as well as some background galaxies and a 
few stars, are clearly evident. 

We performed object detection using SExtractor (Bertin 
& Arnouts 1996) with an RMS weight image that included 
the photon noise and SBF contributions from the subtracted 
galaxy. The full procedure is described in detail by Jordan 
et al. (2004, 2007) and Barber DeGraaff et al. (2007). For the 
most part, we use F814W to define the sample limits because 
it is a broad bandpass at the red end of the optical spectrum, 
and the signal-to-noise (S/N) of the data is high, despite the 
half-orbit integration. For object detection in the ACS im- 
ages, we required an area of at least four connected pixels 
above a S/N threshold of two; thus, a minimum S/N of four 
within the isophotal detection area. The SBF is quite strong in 
the WFC3 F160W image, and we therefore set the detection 
threshold higher, requiring a total S/N of at least 6 to avoid 
spurious detections. Objects were detected and measured in 
each band separately (rather than "dual image mode"), then 
the catalogues were matched across bands; this process re- 
moves spurious objects from the combined multi-band sam- 
ples. Figure |2] illustrates the depth of the detection in F814W 
by plotting the magnitude error within the isophotal area as 
a function of the SExtractor MAG_AUTO parameter, the esti- 
mated total magnitude from SExtractor. 

Candidate GCs in NGC 1399 must be nearly point-like ob- 
jects in the expected magnitude range. Villegas et al. (2010) 
found that the turnover in the GC luminosity function (GCLF) 
occurs at z&su = 22.802 ± 0.044, where Z&50 is the total magni- 
tude in F850LP from the ACSFCS. We find in Section that 
the mean color of the GC candidates is /si4-Z850 ~ 0.3, where 




/ 814 (AB mag) 

Fig. 2. — The error in magnitude for pixels above the isophotal de- 
tection threshold is plotted against hi4 MAG_AUTO values from SEx- 
tractor. The magnitudes are on the AB system. The vertical dashed 
line at /sm = 23.5 mag is drawn somewhat fainter than the mean mag- 
nitude (or "turnover") of the GC luminosity function, which occurs 
at 7gi4 = 23.5 mag (see text). The median isophotal detection error at 
7 8 i4 = 23.5 is 0.028 mag. 
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Fig. 3. — Full-width half maximum (FWHM) is plotted against /gu 
MAG_AUTO parameter from SExtractor for all objects detected in the 
ACS/WFC F814W image. The FWHM values are in pixels, at a scale 
of 0."05 pix~' . The gray shaded region shows the initial selection for 
GC candidates in this plane. The turnover of the GCLF occurs near 
/8i4 ~ 23.2 mag. 

7gi4 is the SExtractor MAG_AUTO value; this would indicate 
a turnover in F814W of ^23.1 mag. This is consistent with 
a magnitude histogram of the GC candidates selected below, 
although based on the Virgo GCLF and the relative distances 
of the two clusters, the GCLF turnover would be expected to 
occur about 0.2 mag fainter. In any case, this is still brighter 
than the dashed line shown at /si4 = 23.5 mag in Figure[2] 

Figure [3] shows the full width at half maximum (FWHM) 
values measured with SExtractor as a function of 7si4 mag- 
nitude. There is a "finger" of compact sources with FWHM 
~ 2 pix that can be distinguished from the background popu- 
lation down to 7gi4 24 mag (shaded region). This is similar 
to the FWHM of the point spread function, which is about 
0"09, or 1.8 to 1.9 pix. At the 20 Mpc distance of Fornax, 
1" corresponds to about 96 pc, so the typical GC half-light 
radius of r/, ss3 pc corresponds to 0.6 pix for the ACS/WFC. 
Thus, the GCs are very marginally resolved, and some with 
smaller r/, will be indistinguishable from point sources. For 
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our initial GC candidate selection, we therefore selected all 
objects in the magnitude range 19.5 < 7si4 < 23.5 mag, with 
1 < FWHM < 4. We also require the candidates to be rea- 
sonably round, with ellipticity < 1/3 (in practice, this rejects 
very few objects, since they are already required to be com- 
pact). In Section [6] where we match our sample against that 
of the ACSFCS, we confirm that the overwhelming majority 
of the objects are high-probability GCs. 

The catalogs produced by SExtractor include magnitudes 
measured within many different apertures. In general, we 
found that the magnitudes of GC candidates measured within 
the SExtractor 6 pix (diameter) aperture proved to be optimal, 
in the sense that the scatter between colors was near the min- 
imum, yet the aperture enclosed ^75% of the total light for 
point-like objects, thus making a good compromise between 
statistical and systematic errors (cf. Appendix F of Sirianni 
et al. 2005). This proved to be the case for both the ACS/WFC 
and WFC3/IR photometry. Although for the ACS/WFC this 
aperture corresponds to a radius r = Of 15, while for the driz- 
zled WFC3/IR data it is r = Of 30, the FWHM of the point 
spread function (PSF) of WFC3/IR is roughly twice that of 
ACS/WFC; thus the aperture corrections from Sirianni et al. 
(2005) and Kalirai et al. (2009) are very similar when the 
WFC3/IR aperture is twice that of the ACS/WFC aperture in 
angular units. As a result, the systematic offsets in the aper- 
ture colors due to differential PSF effects across bands are 
relatively small. 

Following Jordan et al. (2009), we calculated the color 
corrections due to differential aperture effects within our 
apertures (r = Of 15 for ACS/WFC; r = 0."30 for WFC3/IR) 
for a typical GC (King model with rj, = 3 pc and con- 
centration c = 1.5) convolved with the PSF. For £475-/814* 
y&xr hiAi V(,o(,-Hh,o, and /8i4-#i60, tne estimated corrections 
are: +0.02, +0.01, -0.01, and -0.02 mag, respectively. Since 
these offsets are small, systematic, and systematically uncer- 
tain at the ~ 0.02 mag level, we have not applied them to 
our colors for this analysis, but simply note that such offsets 
would be expected for some external comparisons. The con- 
clusions of this work would not change. In contrast, the cor- 
responding corrections would be ~ 0.1 mag for colors involv- 
ing Z850, since the PSF is significantly broader in that band; 
however, in that case (see Section |6j, we use the ACSFCS 
photometry, for which the colors have been converted to the 
infinite-aperture values assuming the King model profile for a 
typical GC convolved with the PSF (see Jordan et al. 2009). 

4. COLOR DISTRIBUTIONS 

NGC 1399 was one of the first galaxies suggested to show 
evidence for bimodality in its GC metal distribution, based on 
Washington system photometry (Ostrov et al. 1993), although 
there was no clear separation between the proposed compo- 
nents in the color distribution of that early sample. The bi- 
modality result, again based on photometric colors, was more 
evident for some other galaxies (e.g., Zepf & Ashman 1993). 
It is now clear that the situation in NGC 1399 is more com- 
plex. Color bimodality is present in the system, but not for all 
subsamples, and we first address this issue. 

4.1. Optical Colors by Magnitude 

Dirsch et al. (2003) discussed the "striking" difference in 
the color distribution of NGC 1399's brightest GCs compared 
to those within < 2 mag of the GCLF turnover, based on wide- 
field CTIO 4 m Mosaic imaging extending to ~ 20'. The GCs 
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FIG. 4. — Optical color-magnitude diagrams for GC candidates in 
NGC 1 399 from ACS/WFC imaging. Dashed horizontal lines indicate the 
range over which the color histograms are significantly bimodal, and the ver- 
tical lines indicate the peaks in the color distributions for this magnitude 
range. The brightest GCs, more than ~2 mag above the turnover, do not 
exhibit distinct bimodality because NGC 1399 is a strong "blue tilt" galaxy. 

in their brightest magnitude bin followed a unimodal color 
distribution peaking at C—T\ « 1.5 in the Washington sys- 
tem, whereas GCs in the next two magnitude bins exhibited 
a clearly bimodal color distribution with the "gap" occurring 
near the same C—T\ ss 1.5 color. This confirmed the earlier 
result based on a smaller sample by Ostrov et al. (1998) that 
the brightest ~ 1 mag of the GC population had a broad color 
distribution, with a mean color similar to the location of the 
gap in the distribution of the fainter GCs. Dirsch et al. also 
found that the red peak was only prominent at r < 9', and is 
more of a red tail at larger radii. Bassino et al. (2006) analyzed 
two additional Mosaic fields to extend the areal coverage of 
NGC 1399 GC system even further, and found results similar 
to those of Dirsch et al. (2003). 

Figure|4]shows 0^^475-/814 and V(,o6~hi4 versus 7s 14 color- 
magnitude diagrams (CMDs) for GC candidates selected in 
7gi4, as described above, and matched with objects detected in 
the GO- 1091 1 £475 data (taken at the same pointing and within 
the same orbit as 7gi4) and, separately, with objects detected 
in V<506 from the multi-pointing GO-10129 observations. The 
colors are measured within the adopted r = 3 pix aperture, and 
we only consider objects with color errors < 0.2 mag for this 
aperture. Many large galaxies exhibit a "blue tilt" in their 
CMDs, a tendency for GCs in the blue peak to become redder 
at higher luminosities, sometimes merging together with the 
red component (Harris et al. 2006; Mieske et al. 2006; S trader 
et al. 2006; Peng et al. 2009). In this sense, NGC 1399 is 
a very strong "blue tilt galaxy," since the two color compo- 
nents merge for the brightest GCs (Dirsch et al. 2003; Mieske 
et al. 2010). However, Forte et al. (2007) find that for fainter 
GCs, in the luminosity range where the distribution is clearly 
bimodal, there is no significant slope in the color of the blue 
peak with magnitude. This is consistent with models in which 
the tilt is due to self-enrichment and requires some minimum 
GC mass threshold (Bailin & Harris 2009). Blakeslee et al. 
(2010a) showed that a simple mass-metallicity relation, com- 
bined with a nonlinear color-metallicity transformation, can 
produce bimodal color distributions with a blue tilt in GC 
populations with broad unimodal metallicity distributions. In 
Figure |U the horizontal lines enclose the magnitude range 
in NGC 1399 where the bimodality is most evident, and the 
vertical lines show the approximate locations of the peaks, 
found below. The tilting of the blue GCs towards the red at 
7814 < 21.5 is clear. 
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Fig. 5. — Histograms of 5475-/814 colors for different ku magnitude ranges. 
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Fig. 6. — Histograms of Veo6-hi4 colors for different /sw magnitude ranges. 
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FIG. 7. — Optical-IR color-magnitude diagrams. The dashed horizontal 
lines are shown to facilitate comparison with the 7gx4 region enclosed by sim- 
ilar lines in Figure|4] they have been shifted by the mean 7gi4— i/i60 color of 
0.70 mag. 

Figures [5] and [6] display #475-/814 and V<j06— isw histograms 
for GC candidates in several different 7si4 magnitude ranges 
(note that there are only four brighter than 7gi4 = 20). In both 
cases, the appearance of bimodality is greatest within the 
21.5 < /gw < 22.5 and 22.5 < I&u < 23.5 magnitude ranges. 
Although the #475-/814 colors in the brightest magnitude range 
of Figure|5]appear somewhat bimodal, the blue "peak" in this 
bright range coincides with the "gap" at #475-/814 ~ 0.95 in 
the color histograms of the fainter GCs. The situation is sim- 
ilar in Figure |6j although in that case, the histogram of the 
brightest GCs is essentially a broad red distribution that en- 
compasses the gap at Veo6~ hu = 0.4 for the fainter GCs. 

Thus, although the Dirsch et al. (2003) and Bassino et al. 
(2006) studies covered a much larger area than our HST data, 
we find the same lack of obvious bimodality in the color dis- 
tribution of the brightest GCs. Our principal goal here is 
to examine the optical-IR colors for GCs that define a dis- 
tinctly bimodal distribution in purely optical colors. This pro- 
vides a simple, fairly direct test of whether the colors linearly 
trace metallicity and reflect true bimodality in the underly- 
ing distribution. In the following section, we therefore com- 
pare the purely optical and optical-IR color distributions over 
the 21.5 < /814 < 23.5 magnitude range where the optical bi- 
modality is strongest. 

4.2. Optical-IR Color Distributions 

For comparison to Figure [4] we show the V6O6-Z/16O and 
/814-//160 versus //160 CMDs in Figure|7] The colors are mea- 
sured within apertures of r=3 pix, as described in Section[3] 
Again, only objects with color errors < 0.2 mag are shown. 
The horizontal lines here are the equivalent of those shown in 
Figure H but shifted by the mean (hu-H m ) = 0.70 ±0.01 
mag. They are merely for comparison to the previous CMDs, 
not to select a magnitude range with the most bimodality. Any 
bimodality in these CMDs is much less evident. 

Figure |8] shows the binned #475-/814, Vm-hu, h\4~H m , 
and V6O6-//160 color distributions, along with smooth density 
estimates constructed with a Gaussian kernel. The colors in- 
volving /814 use the magnitude range selected above, while 
for V606 - Z/i60, we use 22 < V606 < 24, as the mean V606-/814 
color is 0.49 mag. We did this for convenience because the 
pairs of bandpasses were matched separately at this stage; in 
Section [6] we examine the colors of a single merged sample. 
In any case, the selection on either /si4 or V606 makes no dif- 
ference to the conclusions here. The #475-/814 and V606-/si4 
color distributions appear strongly bimodal, both having two 



distinct peaks with the red peak containing roughly 2/3 of the 
GC population. The predominance of red GCs occurs be- 
cause the HST fields cover a relatively small central area, and 
the NGC 1399 GC system has a radial color gradient. The 
/814-//160 and V6O6-Z/16O distributions appear quite similar to 
each other: unlike in the optical colors, there are no well- 
separated peaks, although neither do they appear to be sim- 
ple Gaussians. The V506-//160 distribution shows some evi- 
dence for a blue component; this is understandable because 
/814-//160 is sharply peaked (though, again, not symmetric), 
while Vooerhu lS bimodal, so Vooer- Hi60 would be expected 
to show some bimodality. 

We have used the GMM code ("Gaussian Mixture Mod- 
eling") of Muratov & Gnedin (2010) to quantify the above 
qualitative impressions. These authors provide a detailed dis- 
cussion of the issues and pitfalls inherent in bimodality tests. 
They note that Gaussian mixture modeling such as GMM or 
KMM (Ashman et al. 1994), by itself, is more a test of Gaus- 
sianity than unimodality. As observed by many authors, uni- 
modal but skewed distributions will strongly favor a double 
Gaussian model, even if there is no true bimodality in the dis- 
tribution. To address this point, they emphasize the impor- 
tance of the kurtosis, stating that "kurt < is a necessary but 
not sufficient condition of bimodality." This is true because 
the sum of two populations with different means is necessar- 
ily broader than a single population; therefore, in the case of 
Gaussians, which have kurt = 0, the kurtosis must be nega- 
tive for a valid double Gaussian decomposition. In addition, 
following Ashman et al. (1994), Muratov & Gnedin define 
the quantity D as the separation between the means of the 
component Gaussians relative to their widths; they state that 
Gaussian splits that appear significant based on the p value 
but have D < 2 are "not meaningful," in the sense that they do 
not demonstrate clear bimodality, but non-Gaussianity. 

Table [2] presents our GMM analysis results; the errors on 
the tabulated quantities come from the bootstrap resampling 
done in the GMM code. Note that the analysis is indepen- 
dent of any binning. The full #475-/814 sample for the 21.5 < 
/814 < 23.5 magnitude range contains 584 objects and strongly 
prefers a double Gaussian model with peaks at 0.83 ±0.02 
and 1.19 ±0.01 mag. (We remind the reader that the small 
offsets for aperture corrections, as given in Section [3] should 
be applied for external comparisons.) About 70% of the GC 
candidates are assigned to the red component. The kurtosis 
is negative, and D is significantly above 2, confirming a valid 
bimodal Gaussian model. We therefore assign a "Y" in the 
last column of Table|2]to indicate the affirmative on the ques- 
tion of bimodality. However, because the #475-/814 distribu- 
tion has a few outliers in tails, and mixture modeling codes 
are generally sensitive to extended tails, we ran another test 
after restricting the range to 0.6 < #475-/814 < 1.5. The re- 
sults, given in the second row of Table [2] are very similar 
to prior run, but with the kurtosis being even more negative 
and the model uncertainties slightly reduced. The situation is 
very similar for Vaoer- hu< which favors a well-separated dou- 
ble Gaussian model with peaks at 0.34 ± 0.01 and 0.54 ± 0.01 
mag, and 74 ± 3% of the GCs assigned to the red peak. Again, 
we answer "Y" to the question of bimodality. 

Consistent with the visual impression, the GMM code finds 
that the /814-//160 distribution is not bimodal. When run on 
the full distribution in Figure [8] the few objects in the tails 
result in a positive kurtosis and virtually all of the objects are 
assigned to the red peak. When we restrict the color range 
to 0.1 < /814-//160 < 1.3, the kurtosis becomes negative as 
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Fig. 8. — Histograms of £475-/814, V606-/si4, /si4-f/i60, and V606-/?i60 colors for all matched GC candidates in the specified magnitude ranges 
and with color errors less than 0.2 mag. 



required, but the p value of the double Gaussian model is not 
significant, and the fraction of objects assigned to the second 
component is consistent with zero. 

For V606~#i60> the situation is a bit more nuanced. For GC 
candidates in the broad color range < Veoe-Z/ieo < 2, the 
kurtosis is negative, and the p value of 0.01 is marginally sig- 
nificant, but the separation D is not meaningful. Moreover, the 
weakly favored decomposition has only 39 ± 22% of the GCs 
in the red peak, as compared to <~ 70% for the optical decom- 
position. If we remove the two bluest outliers, the p value be- 
comes slightly less significant, but now D is marginally above 
2, and the fraction of objects in the second (red) peak is con- 
sistent with that found in the optical, although it is also con- 
sistent with 100% at the 2-cr level. We therefore answer an 
uncertain "Y?" to the question of bimodality here, consistent 
with unimodality in /814-//16O but bimodality in V606-4i4- 

As another test of the bimodality in the color distributions, 
we ran the "Dip" test (Hartigan & Hartigan 1985) supplied 
with GMM, which nonparametrically measures the signifi- 
cance of any gaps, or "dips," in a distribution. Since the Dip 
test is insensitive to the assumption of Gaussianity, it is much 
more robust, but as discussed by Muratov & Gnedin, the re- 
turned significance levels are much lower than for the para- 
metric GMM analysis. Probabilities above <~ 50% may be 
taken as indicative of likely bimodality. From this test, we 
found that the probability for the most significant gaps in the 
#475-^814, V606-/814, k\4-H m , and V6O6-//16O distributions to 



be real were 75%, 97%, 4%, and 22%. 

Overall, we confirm that the £475-/814 and Veoe-Zsu distri- 
butions for GCs in NGC 1399 are strongly bimodal, with con- 
sistent proportions of red and blue GCs, when the brightest 
GCs are excluded. The /814-//16O distribution is not signifi- 
cantly bimodal, but there is marginal evidence for bimodal- 
ity in V606-H)60 consistent with the proportions found for 
#475-/814 and V606 _ /8i4- This is understandable if the opti- 
cal color bimodality is partly the result of a changing hori- 
zontal branch morphology, since the colors involving bluer 
bandpasses would be more affected by the behavior of the 
hot horizontal branch stars, whereas colors such as /814-//160 
would not. Thus, some evidence for bimodality in V6O6-//16O 
would be expected, but weaker than that found for V606-/314 
(Cantiello & Blakeslee 2007). 

5. COLOR-COLOR CURVATURE 

The broadest baseline color-color combination we have 
available for investigating stellar population issues is #475-/814 
versus /8i4-£/i60- These color indices span factors of 1.7-2.0 
in wavelength. The #475-/814 index for GCs is affected by 
the properties of stars near the main sequence turnoff, on the 
horizontal branch, and on the giant branch. It is therefore 
sensitive to metallicity, age, and any other parameters (e.g., 
helium) that affect the temperature of the main sequence and 
the morphology of the horizontal branch. By comparison, the 
V606-/814 index spans a much smaller factor (1.3) in wave- 
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Fig. 9. — The ACS+WFC3/IR hu—Hieo versus #475-/814 color-color plots. In the left panel, the invidual matched GC candidates are shown, 
with error bars estimated from the RMS images mentioned in Section [3] The red curve represents a 4th-order robust polynomial fit with 
coefficients given in the text. The right panel shows the median values within ten bins (40 objects per bin) ordered by #475-/814. We also plot 
the same polynomial fit as in the left panel (red dashed curve), a 4th-order polynomial fit to the median points (blue curve), and a linear fit to 
these points (black dashed line). Error bars on the points are plotted, but in most cases they are smaller than the point size. The plot scale is 
expanded in the right panel to clarify the differences between the fits. 



length, and is less sensitive than #475-/814 to features such as 
the blue horizontal branch. As discussed in the Introduction, 
colors such as /8i4~#i60 are mainly sensitive to the temper- 
ature of the giant branch, which is controlled by metallicity 
with very little age dependence. Colors such as gAi5~H\a) 
and V606 - #160 are aiso sensitive to giant branch temperature, 
but they complicate the metallicity dependence with addi- 
tional sensitivity to the main sequence turnoff and horizontal 
branch. We therefore concentrate now on the relation between 

g415~hl4 and/8i4-//i60- 

Figure [9] (left panel) plots /814— ffi6o as a function of 
g475-/gi4 for 401 NGC 1399 GC candidates in our combined 
ACS-WFC3 data set with 19.5 < / 814 < 23.5 and 0.5 < 
^475-^814 < 1-6. The corresponding data are listed in Table[3] 
Although we deal throughout this work with aperture colors 
within a 3-pixel radius, in order to facilitate external compar- 
isons, the last two columns of Table[3]provide colors that have 
been corrected for differential aperture effects, assuming the 
corrections for a typical GC, as given in Section [3] Since we 
are interested here in the relationship between these color in- 
dices, rather than simply the presence or absence of bimodal- 
ity, we improve the statistics by including the GCs from the 
brightest magnitude range in Figure [5] increasing the sample 
by 17%. The relation appears nonlinear. The plotted curve 
represents a quartic (4th order) polynomial, which has been 
fitted to the data using robust orthogonal regression (Jefferys 
et al. 1988). This approach minimizes the residuals in the di- 
rection orthogonal to the fitted relation, and it is appropriate 
here because of the significant observational scatter in both 
coordinates. The fitted coefficients are given by 

hu-Hm = -14.54 + 49.23x - 59.02x 2 + 30.14x 3 

-5.24x 4 , (1) 

where x = (gxis—hu)- Several of the objects appear discrepant 
with respect to this relation: these may be due to contami- 
nation in the sample from Galactic stars and/or background 



galaxies, or to real stellar population variations. In particu- 
lar, objects that lie above or to the left of the relation could be 
young clusters, or they may have extreme horizontal branches. 
Despite some outliers, there is a tight locus of objects that de- 
fine the curved sequence in color-color space traced by the 
above fit. 

To further illustrate the curvature, we have binned the data 
by ^475-^814 into 10 groups of 40 objects, and we plot the 
median values of the colors within the bins in the right panel 
of Figure [9] The polynomial fit from the left panel is plotted 
again, as well as linear and quartic fits to the binned data. The 
deviation from a linear relation is highly significant. If we use 
the scatter to estimate the errors (a/y/N) in the medians, the 
value of the reduced xl is 2.7 for the linear fit, and 0.1 for 
the quartic fit. Thus, the linear fit is strongly rejected. The 
low value of xl for the quartic fit probably indicates that the 
uncertainties in the medians have been overestimated, most 
likely because the small fraction of outliers make the scatter 
estimate too high. This also suggests that xl should be larger 
than 2.7 for the linear fit. The 4th order model closely traces 
the curvature of the binned points; adding another term to the 
polynomial does not reduce xl further. The polynomial fits 
to the binned and unbinned data diverge near the endpoints 
because of the lack of constraints, but they agree well over 
the ~ 0.8—1 .3 color range. 

Figure [lO] shows the fit residuals as a function of the mean 
of the #475-/814 and h\\-H\(& colors, which corresponds to 
half of the £475— ffi60 color. Since the color-color relation 
has an average slope of about one (the linear fit in Figure [9] 
has slope 1.01 ±0.06), this mean color is, to a good approx- 
imation, proportional to the distance along the relation (and 
therefore in the direction orthogonal to the residuals). The 
small squares in this figure represent individual GC candi- 
dates, and their residuals are with respect to the fit given in 
equation (UJ). Open squares are used for the faintest magnitude 
of GC candidates included in the fitting, while solid squares 
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FIG. 10. — Color-color fit residuals are plotted as a function of 0.5 X 
(§475 _ ^l6o). which corresponds to the mean of £475-/814 and /gi4— //i60- We 
use this mean for the horizontal axis because it is approximately orthogonal 
to the residuals that were minimized in the robust regression fitting (see text). 
Square symbols show the orthogonal residuals for individual GC candidates 
with respect to the curve in the left panel of Fig. [9] Filled squares are used 
for objects with /gi4 < 22.5, while open squares are used for objects with 
22.5 < 414 < 23.5. The large circles indicate the quartic fit residuals for the 
median points shown in the right panel of Fig. [9] 

are used for brighter objects. The fainter objects scatter more, 
but otherwise do not differ systematically from the brighter 
ones. The large circles in Figure QJj] show the residuals for 
the median-filtered points with respect to the quartic fit for 
these points plotted as a solid curve in the right panel of Fig- 
ure [9] In this case, the median values have very little scatter, 
and the fit is essentially identical regardless of whether or- 
thogonal regression or a simple least-squares approach (min- 
imizing residuals in hu-Hi^) is used. From Figure [10] we 
conclude that the form of the color-color relation is well ap- 
proximated by a quartic polynomial over a range in #475-//i60 
from ~ 1 .3 to ~ 2.2 mag, which corresponds approximately to 
the 0.8 < #475-/814 ^1-3 mag range quoted above. 

We emphasize that not only is the relation between 
#475-/814 and /814— Z/16O nonlinear, but the slope attains a local 
minimum within the color range where it is particularly well 
constrained. In fact, the inflection point of the quartic fit to 
the unbinned data occurs at #475-/814 = 1.002, which corre- 
sponds closely to the dip in the #475-/814 color distribution in 
Figure [8] If /814-//160 is a good, relatively simple, indicator 
of metallicity, in line with theoretical expectations, then this 
type of "wavy" or inflected relation could produce double- 
peaked optical color histograms from metallicity distributions 
that have a very different form. The metallicity distribution 
does not need to be unimodal or even symmetric; chemical en- 
richment scenarios tend to produce asymmetric distributions 
(see Yoon et al. 201 lb for a detailed discussion of this issue). 
We address the metallicities briefly in the following section. 

The observed wavy relation between #475-/814 and 
/814-//160 accounts for the marked difference between the 
purely optical and optical-IR color distributions found in Fig- 
ure [8] The significance of the curvature demonstrates that the 
apparent disparity was not the result of underestimated pho- 
tometric errors obscuring bimodality in /8i4 _ //i60- 

6. ACSFCS COMPARISON: IMPLICATIONS FOR METALLICITY 

The ACSFCS (Jordan et al. 2007) is a survey with ACS in 
the F475W and F850LP bands of 43 early-type Fornax clus- 
ter galaxies. All objects with colors and magnitudes con- 



sistent with GCs in the program galaxies were fitted with 
King (1966) models using the methodology of Jordan et al. 
(2005) to derive half-light radii r/, and total magnitudes. As 
described in Jordan et al. (2009), objects with colors in the 
range 0.5 < #475-2850 < 1 -9 and radii 0.75 < r/, < 10 pc were 
assigned probabilities pec of being GCs based on their zsso 
magnitudes and half-light radii, as compared with those found 
for objects in background fields. Those with pec > 0.5 were 
accepted as likely GCs. For the #475-2850 colors considered 
in this section, we used the NGC 1399 photometric catalogue 
produced as described in Jordan et al. (2009) and already used 
in several ACSFCS publications (Masters et al. 2010; Mieske 
et al. 2010; Villegas et al. 2010; Liu et al. 201 1). 

We used the object positions to match the ACSFCS pho- 
tometry for NGC 1399 with our #475, /814, and H160 mea- 
surements. We keep high probability GCs with pec > 0.9 
(in practice, there was only one matched object with 0.5 < 
Pec < 0.9). This merged data set provides two important ben- 
efits. First, the probabilistic selection based on r/, and 2sso 
should remove most of the remaining contaminants in our 
sample. Second, the ACS #475-2850 colors have been cali- 
brated empirically against metallicity. Although the through- 
put for objects with GC spectra is about twice as high in 
F814W as in F850LP, allowing more precise color measure- 
ments with /814 for a given exposure time, the longer base- 
line afforded by zsso improves the metallicity sensitivity (Cote 
et al. 2004). Peng et al. (2006) found that the relation be- 
tween [Fe/H] and #475-2850 was not adequately described by 
a linear fit; they used a broken linear model with a shal- 
lower slope for #475-2850 > 105. This was adequate over 
the range #475-2850 = 0.7-1 .4, but under-predicted the spectro- 
scopic metallicities of M87 and M49 GCs in the 1.4-1.6 color 
range. Additional curvature was needed to match these. To 
accommodate the measured metallicities of these redder GCs, 
Blakeslee et al. (2010a) used a polynomial fit that followed the 
data over the full 0.7-1.6 range in #475-2850- This empirical 
calibration provides us with some handle on the metallicities 
for our sample. 

Figure [TTJ shows the color and estimated metallicity his- 
tograms for our sample of ACS+WFC3 data after merging 
with the ACSFCS dataset. Again we use the 21.5 < 7gi4 < 
23.5 magnitude range for the histograms to select the regime 
of optical color bimodality. After the merging with ACS- 
FCS, there was one object at #475-/814 = 0.45 mag, whereas 
all the others were in the 0.64-1.52 range; we removed this 
single outlier in #475-/814 color. In contrast to the approach in 
Figure [8] the exact same 322 high-probability GCs are repre- 
sented in all of the color histograms in Figure [TTJ For unifor- 
mity, we have used the same binsize of 0.07 mag for all the 
colors. The change in binning is mainly responsible for the 
different appearance of the Zsw-Z/ieo histograms in Figures[8] 
andHU but the kernel density curves reflect the same features. 

As described above, the metallicity histogram in the lower 
right of Figure[TT|is based on the polynomial fit by Blakeslee 
et al. (2010a) to the Peng et al. (2006) data. We use this poly- 
nomial transformation only for objects with #475-2850 < 1-6 
mag, since it is not constrained beyond this; objects with red- 
der #475-2850 colors (7.5% of our merged sample) are likely 
affected by observational scatter. Thus, there are fewer ob- 
jects in the [Fe/H] histogram, and the truncation at [Fe/H] = 
0.76 dex is artificial. However, the peak at [Fe/H] ss -0.3 dex, 
and the long tail to lower metallicities, accurately reflect the 
#475-2850 distribution combined with the empirical metallicity 
calibration. 
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Fig. 11. — Histograms of #475-/814, £475-^850, and hu-Hi6o colors, along with [Fe/H] values inferred from ACSFCS g475-Z850 data (see text), 
are plotted for GC candidates identified in the ACS Fornax Cluster Survey. 



Table|4]summarizes the GMM results for this merged sam- 
ple. For the #475-/814 distribution, we again find a signifi- 
cant preference for a double Gaussian model, with peaks at 
0.84 ±0.02 and 1.20 ±0.01 mag and 73 ± 5% of the objects 
in the red component. The result is similar for £475— Z850, 
with peaks at 0.98 ± 0.04 and 1.41 ±0.01 mag and 72 ± 7% 
in the red component. When the analysis is run on the full 
/gi4-ffi6o sample, we find kurt > 0, D < 2, and p > 0.05; 
thus, no significant evidence for a bimodal model. How- 
ever, no outliers have been trimmed from the /8i4~//i60 dis- 
tribution; if we remove the two objects with hu-H\eo < 0, 
then the fourth row of Table [4] shows kurt < 0, as required 
for significant Gaussian bimodality, but the p value increases 
to 0.055, and D is unchanged, except in the uncertainty. If 
we trim the outliers more aggressively, restricting the range 
to 0.2 < hu-Hi(,o < 1.2, then kurt decreases further, and 
the separation between the best-fit peaks is now a marginally 
significant D = 2.0 ±0.3, but p increases to 0.062 (signifi- 
cance less than 2a). Moreover, the preferred breakdown has 
24 ± 28% of the objects in the red component (before the trim- 
ming, the red percentage was even lower). This red fraction is 
consistent with zero; however, if it were interpreted as real, 
then this would indicate a "different bimodality" than that 
found for the optical colors, where ^73% of the objects are 
in the red peak. Thus, the preferred optical and optical-IR 
bimodalities could not both reflect a common, more funda- 
mental bimodality in metallicity. 



We also experimented with the GMM analysis for the 
[Fe/H] distribution derived from the £475-2850 colors. For the 
full sample, the kurtosis is significantly above zero, and the 
code wants to fit a primary peak with 80% of the objects, 
then use another broad Gaussian for the blue tail. This re- 
flects the obvious lack of Gaussianity. If we aggressively trim 
the blue tail, removing objects with [Fe/H] < -1.9, then we 
find kurt < 0, p = 0.002, and D > 2. However, the fraction 
of objects in the main peak is then 0.917 ± 0.084, consistent 
with 100%. More to the point, this procedure of removing 
members of a continuous extended distribution (rather than 
outliers), is not statistically valid; we have done so only for il- 
lustration. The metallicities for NGC 1399 GCs derived from 
the (bimodal) £475-2850 colors and empirical calibration de- 
fine a continuous distribution with a peak at ~ -0.3 dex and 
an extended blue tail. The form of the high-metallicity end is 
unknown, but it is unlikely to extend much beyond the limit 
imposed here by the calibration. 

Figure[l2]plots ku-H m versus £475-/814 and g 4 75-Z850 for 
the merged sample. This includes 385 objects with 19.5 < 
7gi4 < 23.5; thus, the merging with ACSFCS removed only 16 
objects, or 4% of the sample. The resulting plot for #475-/814 
is essentially identical to that in Figure|9] The robust orthogo- 
nal regression insured that the rejected objects (possible con- 
taminants) did not strongly affect the fit. To illustrate the con- 
sistency between the two relations in Figure [12] we have not 
independently fitted the relation for the ACSFCS £475-2850 
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Fig. 12. — Similar to the left panel of Figure [5] but here shown for the merged sample of high-probability GC candidates identified in the 
ACSFCS. The g475-Z850 measurements are from the ACSFCS. 



colors, but simply transformed the fit for #475-/814 aperture 
colors according to the linear relation we find for the data: 

(#475-^8 14)r=3p = 0.13 + 0.75(£ 475 -285o)aCSFCS • (2) 

It is not surprising that the relations between /8i4~//i60 an d 
these two optical colors would be very similar, since they have 
the same blue bandpass, and there is substantial wavelength 
overlap of the red bandpasses. We note however that the #475 
measurements are from independent data sets, one taken at the 
same time and orientation as the zsso data in HST Cycle 13, 
and the other taken with the Igu data in Cycle 15. 

The merged sample analyzed here is nearly 5 times larger 
than previous optical-IR color studies of GCs with HST, and 
has significantly better precision than previous ground-based 
studies. Clearly, even larger, high-quality optical-IR samples 
are needed to explore higher-order features of the color-color 
relations presented here. Deep spectroscopy of representa- 
tive portions of the photometric samples is also essential for 
more detailed constraints on the underlying metallicity dis- 
tributions. We are currently analyzing the full sample of 16 
galaxies from our WFC3/IR program, all of which have exist- 
ing ACS £475-2850 data. The results from that larger sample 
will allow us to test the conclusions from the present study. 

7. DISCUSSION AND CONCLUSIONS 

In this study, we have combined ACS and WFC3/IR pho- 
tometry to examine optical and optical-IR color distributions 
and color-color relations for —400 GCs in NGC 1399, the 
central giant elliptical in Fornax. Consistent with previous 
wide-field studies (e.g., Dirsch et al. 2003), we find that 
the brightest GCs (about 2-3 mag brighter than the GCLF 
turnover) in NGC 1 399 have a broad red distribution in optical 
colors. Fainter than this, the optical colors become strongly 
bimodal, with — 70% of the GCs in the red component. Any 
bimodality in the optical-IR colors is much less clear; for 
hi4~H\(,o, the preferred relative proportion of red and blue 
GCs also differs from in the optical. 

In addition, the relation between pure optical and optical- 
IR color is significantly nonlinear. The slope of the depen- 
dence of /814-//16O on £475-/814 has a local minimum near 



#475-/8 14 ~ 1 (the polynomial fit has a true inflection point 
here), corresponding approximately to the dip in the #475-/814 
distribution; the result is similar for £475-2850- Since GC 
color is mainly determined by metallicity, at least one of the 
color-metallicity relations must also be significantly nonlin- 
ear. Such nonlinearities can produce bimodal color distribu- 
tions even when the metallicity distribution lacks any clear 
bimodality (Richtler 2006; Yoon et al. 2006; Cantiello & 
Blakeslee 2007). The main requirement is that the slope in 
metallicity as a function of color must have a local minimum 
near the mid-range of the color distribution (for additional dis- 
cussion, see Blakeslee et al. 2010a). Such "wavy" nonlinear- 
ity for optical colors is predicted by some stellar population 
models (Lee et al. 2002; Yoon et al. 2006, 201 la,b; Cantiello 
& Blakeslee 2007), largely as a result of the behavior of the 
horizontal branch (still imperfectly modeled). 

Further evidence for nonlinearity comes from the observed 
relation of £475-2850 colors of GCs with respect to spectro- 
scopically estimated metallicities (Peng et al. 2006). We have 
used a smooth fit to this relation, along with £475-2850 col- 
ors from the ACSFCS, to estimate metallicities for the GCs in 
our sample. Matching with the ACSFCS, which used an in- 
dependent, probabilistic method for selecting GCs, also con- 
firms that the level of contamination in our sample is quite 
small. The metallicities derived from the £475-2850 colors fol- 
low a continuous distribution that peaks at — -0.3 dex and 
extends in a low-metallicity tail to [Fe/H] w -3. Because 
of the limited range of the calibration, the form of the high- 
metallicity end remains uncertain, but it cannot have the same 
long extent in logarithmic space as the low-metallicity tail. 
This skewed distribution, based on the empirical calibration, 
closely resembles the GC metallicity distributions derived by 
Yoon et al. (2011b) using their model color-metallicity rela- 
tions. As these authors have pointed out, the resulting distri- 
butions are in much better agreement with those of the halo 
stars in the same early-type galaxies. 

The relative importance of the form of the color-metallicity 
relation on the one hand, and structure in the underlying 
metallicity distribution on the other hand, remains an unre- 
solved question, at least in the minds of many. One argument 
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against the importance of nonlinear color-metallicity relations 
in extragalactic GC populations comes from the bimodal GC 
metallicity distribution in the Milky Way Galaxy (e.g., Zinn 
1985). However, it is unclear if the GC system of the Milky 
Way would be recognized as bimodal to an outside observer, 
given the necessity of avoiding disk contamination when se- 
lecting GCs in spirals. For instance, Zepf & Ashman (1993) 
show that the metallicity distribution of the Milky Way halo 
GCs, the subpopulation that they considered relevant for com- 
parison with elliptical galaxies, is unimodal. Furthermore, the 
generalization of the Milky Way's properties to galaxies such 
as NGC 1399, which has ^40 times as many GCs and must 
have had a much more chaotic interaction history than a spiral 
galaxy in a low-density group, is uncertain. It is clear that the 
metallicity peaks in the Milky Way do not correspond with 
those derived for giant ellipticals based on linear transforma- 
tions of optical colors (e.g., Zepf & Ashman 1993; Gebhardt 
& Kissler-Patig 1999; Peng et al. 2006). 

Several other large, nearby galaxies possess extensive GC 
systems that are probably better analogues to those of giant 
ellipticals. The nearest of these is M3 1 , with roughly 500 con- 
firmed GCs (Huxor et al. 201 1). Based on high-quality spec- 
troscopic metallicities for 322 M3 1 GCs (twice as many as in 
the whole of the Milky Way), Caldwell et al. (201 1) find that 
the metallicity distribution does not evince clear bimodality, 
"in strong distinction with the bimodal Galactic globular dis- 
tribution." Instead, M31 shows a broad peak at [Fe/H] w — 1, 
with possibly several minor peaks superposed like ripples on 
the broader distribution. NGC 5128 (Centaurus A) is another 
galaxy with an extensive GC system having a broad, com- 
plex spectroscopic GC metallicity distribution with perhaps 
three closely spaced peaks (Beasley et al. 2008; Woodley et al. 
2010; see discussion in Yoon et al. 201 lb). The evidence sug- 
gests that both of these large nearby galaxies have had more 
active formation histories than has the Milky Way, although 
given their locations in small groups, it is unlikely that their 
histories were as turbulent as those of cluster ellipticals. 

The recent innovative studies by Foster et al. (2010, 2011) 
find that metallicity distributions of GCs inferred from the cal- 
cium triplet (CaT) index differ from those found by simple lin- 
ear transformation of the optical colors. Foster et al. (2011) 
report that the peak in metallicity from their CaT measure- 
ments in the elliptical galaxy NGC 4494 actually corresponds 
to the trough in the g—i distribution, and they note the similar- 
ity of this result to their CaT study of the elliptical NGC 1407 
(Foster et al. 2010), as well as the results on M31 from Lick 
indices by Caldwell et al. (2011). The authors conclude that 
the lack of bimodality in the CaT distributions, and the sim- 
ilarity in the CaT-derived metallicities for bright GCs on op- 
posite sides of the color trough, "casts serious doubt on the 
reliability of the CaT as a metallicity indicator" (Foster et al. 
201 1). However, they also acknowledge that the very different 
color and metallicity distributions in these elliptical galaxies 
can be reconciled if the optical color bimodality is mainly the 
result of nonlinearity in the color-metallicity relation. In this 
case, the trough in the color distribution would correspond to 
a minimum in the slope of metallicity as a function of color; 
thus, GCs straddling the color divide would have nearly the 
same metallicities, as observed. 

Further evidence on both sides of the bimodality debate 
from spectroscopic and optical/near-IR observations of extra- 
galactic GC systems has been discussed by Blakeslee et al. 
(2010a; see Sec. 4 of that work) and in even greater detail 
by Yoon et al. (2011b; see their Sec. 2). We refer the reader 



to those works for more information. The main conclusions 
were that while evidence exists for breaks in the metallicity 
distributions of some galaxies, clear bimodality like that seen 
in the Milky Way is rare. In particular, the distribution in 
metallicity usually appears less bimodal than that of the opti- 
cal colors, which suggests that at least some of the observed 
color bimodality is likely the result of the color-metallicity 
relation. Pipino et al. (2007) arrived at a similar conclu- 
sion, based mainly on the spectroscopic results of Puzia et al. 
(2005). More recently, Alves-Brito et al. (2011) have found 
spectroscopic evidence for bimodality in the metallicities of 
GCs in the Sombrero galaxy (NGC 4594), a nearly edge-on 
Sa galaxy with a prominent bulge. Thus, in contrast to M31 
and Cen A, this may be an example of another galaxy, in ad- 
dition to the Milky Way, with a truly bimodal GC metallicity 
distribution. 

We now summarize the conclusions from this optical-IR 
study of GCs in NGC 1399: 

• The optical g— z, g—I, V-I color distributions are bi- 
modal; the evidence is much weaker for I-H. If one 
is determined to enforce bimodality via double Gaus- 
sian modeling, the resulting bimodal split (dominated 
by the blue component) is different from that found for 
the optical colors (dominated by the red). Thus, these 
color bimodalities cannot both be linear reflections of 
an underlying metallicity bimodality. 

• These "differing bimodalities" imply that there must 
be a nonlinear relation between the purely optical and 
mixed optical-IR colors; indeed, we find empirically 
that the dependence of optical-IR I-H color on optical 
g—I and g— z is nonlinear. The significance of the curva- 
ture indicates that it is not simply a case of photometric 
errors obscuring bimodality in I-H. 

• At least some of the colors studied here must vary non- 
linearly with metallicity, with the slope of the metal- 
licity dependence having a local extremum somewhere 
within the observed range of these data. Observations 
and modeling both indicate that, for the case of GC pop- 
ulations, broad baseline optical colors are more likely to 
exhibit this type of nonlinearity. 

• Thus, the "universal" bimodality in the optical GC 
color distributions of elliptical galaxies is at least partly 
the result of the inflected form of the nonlinear color- 
metallicity relations. 

Determining the full importance of the nonlinear behavior 
of GC colors, and the precise form of the underlying metallic- 
ity distributions, will require further effort in amassing large 
samples of high-quality optical/near-IR photometry and spec- 
troscopy, as well as in improved theoretical modeling. As 
we have done here, it is essential to compare the propor- 
tions in any proposed bimodal decompositions to check for 
consistency among techniques, or across wavebands, rather 
than simply quoting the significance at which single-Gaussian 
models can be rejected. We are currently analyzing the 
optical-IR GC color distributions for the full sample of 16 
galaxies from our HST WFC3/IR project (H. Cho et al., in 
preparation); this and other ongoing projects will continue to 
push back the obscuring veil over extragalactic GC metallicity 
distributions and permit more useful constraints on the forma- 
tion histories of massive early-type galaxies. 
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TABLE 1 

Summary of Observing Programs for NGC 1399 



Instr. & 




Prog 




Dataset 


Bandpass 


Exp. time 




mi a mag b 


Channel 




(GO) 








(sec) 




(AB) 


WFC3/IR 




11712 




IB1H08020 


F160W 


1197 




25.960 H l60 


ACS/WFC 




10911 




J9P305010 


F814W 


1224 




25.937 7 814 


ACSAVFC 




10911 




J9P305020 


F475W 


680 




26.068 g 475 


ACS/WFC 




10129 




J8ZQ0[1-4]010 


F606W 


2108 




26.486 V 606 


ACSAVFC 




10217 




J90X02010 


F850LP 


1220 




24.862 z 85 o 


ACSAVFC 




10217 




J90X02020 


F475W 


760 




26.068 ?475 


a Photometric zero point: the AB 


magnitude for corresponding to one count per second. 








b The symbol used here to denote the AB-calibrated magnitudes and other quantities in the given bandpass. 
















TABLE 2 
















GMM Results for Matched Samples 






Quant 


N 


min 


max 


kurt p-val 


Pi 


Pi 


D 


frac(2) bi? 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) (10) 


g~I 


584 


0.36 


1.60 


-0.58 < 0.001 


0.83 ± 0.02 


1.19 ±0.01 


3.4 ± 0.3 


0.68 ± 0.04 Y 


8-1 


579 


0.64 


1.47 


-0.93 < 0.001 


0.82 ± 0.01 


1.18 ±0.01 


3.8 ±0.2 


0.71 ±0.03 Y 


V-I 


578 


0.25 


0.76 


-0.86 < 0.001 


0.34 ± 0.01 


0.54 ± 0.01 


3.3 ± 0.2 


0.74 ± 0.03 Y 


I-H 


345 


-0.22 


1.36 


0.28 0.266 


-0.22 ± 0.36 


0.69 ±0.10 


5.5 ± 1.8 


1.00 ±0.37 N 


I-H 


342 


0.10 


1.25 


-0.43 0.110 


0.67 ±0.13 


0.94 ± 0.09 


1.7 ± 0.4 


0.08 ± 0.30 N 


V-H 


341 


0.09 


1.86 


-0.25 0.010 


1.01 ±0.13 


1.39 ±0.10 


1.6 ±0.4 


0.39 ± 0.22 N 


V-H 


339 


0.27 


1.86 


-0.54 0.012 


0.88 ±0.18 


1.34 ±0.12 


2.1 ±0.5 


0.62 ± 0.29 Y? 


Note. - 


- Columns list: (1) quantity analyzed; (2) number of objects 


in the analyzed sample; 


(3) minimum and maximum values of the given quantity for the 


GCs in the sample; (4) kurtosis of the distribution; (5) GMM "p value," indicating the significance of the preference 


for a double Gaussian over a single Gaussian 


model (lower p values are more significant); (6) mean and uncertainty of the first peak in the double Gaussian model; (7) mean and uncertainty of the second 


peak in the double Gaussian model; (8) separation D of the peaks in units of the Gaussian a; 


(9) fraction of GC candidates assigned to the second Gaussian 


component (or to the "red peak"); (10) assessment of the evidence for bimodality. 


















TABLE 3 














Matched WFC3/IR and ACS Photometry for GC Candidates 




ID 


RA. 




Dec. 


hu 


(g475-hu)r=3 


(/814-#|6o)/-=3 


(#475 


-7814)cor (/814-#lH))cor 




(J2000) 




(J2000) 


(mag) 


(mag) 


(mag) 


(mag) (mag) 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 




(7) (8) 


14 


54.603067 




35.467057 


22.892 ± 0.023 


1.125 ±0.067 


0.642 ± 0.047 


1.145 0.622 


26 


54.603686 




35.467382 


23.390 ± 0.029 


1.283 ±0.077 


0.831 ± 0.062 


1.303 0.811 


40 


54.603018 




35.464966 


21.741 ±0.012 


1.212 ±0.037 


0.881 ±0.020 


1.232 0.861 


198 


54.605421 




35.464718 


20.974 ± 0.013 


1.277 ±0.040 


0.878 ±0.018 


1.297 0.858 


269 


54.605974 




35.464424 


21.457 ±0.012 


0.924 ± 0.033 


0.473 ± 0.021 


0.944 0.453 


300 


54.608004 




35.467911 


23.063 ± 0.021 


1.198 ±0.054 


0.607 ± 0.041 


1.218 0.587 


346 


54.606921 




35.464733 


22.779 ± 0.019 


1.209 ±0.050 


0.719 ± 0.034 


1.229 0.699 


355 


54.608060 




35.466963 


22.781 ±0.018 


1.220 ±0.049 


0.838 ±0.031 


1.240 0.818 


368 


54.605945 




35.461861 


22.172 ±0.016 


0.891 ±0.039 


0.491 ± 0.030 


0.911 0.471 


Note. - 


- Columns list: ( 1 ) ID number; (2) 


right ascension; (3) declination; (4) SExtractor MAG_AUTO value from F814W image; (5) §475— /gi4 color for r=3 


pix measurement aperture; (6) /8i4-77i60 color for r=3 pix measurement aperture; (7) g475-/si4 


color, with estimated correction for differential aperture effects; 


(8) Igu-H 160 color, with estimated correction for differential aperture 


effects. All magnitudes and colors are corrected for Galactic extinction. Note that the 


analysis and fits presented in this paper are based on the aperture colors 


given in columns (5) and (6). 






(This table 


is available in its entirety in a machine-readable form in the online journal. A portior 


is shown here for guidance regarding its form and content.) 












TABLE 4 
















GMM Results for Merged GC Sample 






Quant 


N 


min: 


max 


kurt p-val 


Pi 


P2 


D 


frac(2) bi? 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) (10) 


8-1 


322 


0.64 


1.52 


-0.71 <0.001 


0.84 ± 0.02 


1.20 ±0.01 


3.6 ±0.3 


0.73 ± 0.05 Y 


g-z 


322 


0.76 


1.79 


-0.90 < 0.001 


0.98 ± 0.04 


1.41 ±0.03 


3.2 ± 0.2 


0.72 ± 0.07 Y 


I-H 


322 


-0.22 


1.25 


0.24 0.052 


0.66 ± 0.30 


0.93 ±0.11 


1.7 ± 1.4 


0.10 ±0.36 N 


I-H 


320 


0.10 


1.25 


-0.45 0.055 


0.66 ±0.12 


0.93 ± 0.08 


1.7 ±0.4 


0.10 ±0.29 N 


I-H 


316 


0.21 


1.18 


-0.65 0.062 


0.62 ±0.11 


0.92 ± 0.08 


2.0 ±0.3 


0.24 ± 0.28 N 


[Fe/H] 


298 


-2.87 


0.76 


0.78 < 0.001 


-1.33 ±0.25 


-0.23 ± 0.06 


2.0 ± 0.6 


0.78 ± 0.08 N 


[Fe/H] 


285 


-1.83 


0.76 


-0.22 0.002 


-1.49 ±0.18 


-0.29 ± 0.07 


3.3 ±0.5 


0.92 ± 0.08 N 



NOTE. — Columns are the same as in Table|2] 



